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ABSTRACT 

We present results of gasdynamics+A^-body constrained cosmological simulations of the Lo- 
cal Supercluster region (LSC; about 30ft. _1 Mpc around the Virgo cluster), which closely mimic 
the real Universe within 100 Mpc by imposing constraints from MARK III catalog of galaxy 
peculiar velocities. The simulations are used to study the properties and possible observational 
signatures of intergalactic medium in the LSC region. We find that, in agreement with pre- 
vious unconstrained simulations, w 30% of the gas in this region is in the warm/hot phase at 
T ~ 10 5 — 10 7 K, and « 40% in the diffuse phase at T < 10 5 K in low-density regions. The X-ray 
emission from the warm/hot gas may represent a small (~ 5 — 10%) but important contribution 
to the X-ray background observed by the ROSAT All-Sky Survey at energies around 1 keV. The 
best prospects for detection of the warm/hot intergalactic medium of the LSC located in fila- 
ments and in the vicinity of virialized regions of groups and clusters are through absorption in 
resonant lines of OVII and OVIII in soft X-rays and in the OVI doublet in UV. If intergalactic gas 
in filaments (p/(p) ~ 1 — 10) is enriched to typical metallicities of > 0.05, the column densities 
of OVI, OVII, and OVIII along a random line of sight near the North Galactic Pole, especially 
near the supergalactic plane, have a significant probability to be in the range detectable by the 
current (FUSE, XMM) and future (Constellation-X) instruments. 

Subject headings: cosmology:theory - large-scale structure - methods: numerical 



1. Introduction 

At the present epoch baryons in stars and 
galaxies constitute only a small fraction of the 
total baryon density predicted by the Big Bang 
nucleosynthesis model (Fukugita, Hogan, & Pee- 
bles 1998). This is in agreement with the Cold 
Dark Matter (CDM) models of structure forma- 
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tion which predict that a large fraction of baryons 
at z w is located in filaments and intergalactic 
medium of groups and clusters, where it is shock- 
heated to temperatures of T ~ 10 5 — 10 7 K (Cen 
& Ostriker 1999b; Dave et al. 1999). 

The detection of the warm/hot intergalactic 
medium (WHIM) and studies of its properties rep- 
resent a challenge. The integrated soft X-ray emis- 
sion of this gas may contribute significantly to the 
observed extragalactic X-ray background (XRB; 
e.g., Croft et al. 2001; Phillips et al. 2001; Voit 
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& Bryan 2001). The best prospects for direct de- 
tection are, however, through absorption or emis- 
sion in the lines of ionic species of heavy elements 
such as oxygen, which requires very sensitive X-ray 
and UV spectroscopy of bright background sources 
(Shapiro & Bahcall 1980; Basko et al. 1981; Ald- 
croft ct al. 1994; Hellsten et al. 1998; Perna & 
Loeb 1998; Fang & Canizares 2000). The recent 
launches of new generation X-ray (Chandra and 
XMM-Newton) and UV (FUSE) satellites, as well 
as access to the UV wavelengths with the Hubble 
Space Telescope (HST), led to first detections of 
absorption in the resonance OVI doublet in UV at 
both low (Tripp & Savage 2000; Tripp et al. 2000) 
and intermediate (Reimers et al. 2001) redshifts. 
The first attempt to detect resonance line absorp- 
tion of heavy metals in soft X-rays with Chandra 
by Fang et al. (2001) yielded only an upper limit, 
but more sensitive studies are currently underway. 

Spectroscopic studies of the warm/hot medium 
are very promising for understanding properties 
and evolution of this important baryonic compo- 
nent. An important piece of information is the 
state of the WHIM gas in the nearby Universe 
at z = 0. However, the absorption at very low 
redshifts should be due to only a few large-scale 
structures and one would therefore expect model 
predictions to be rather uncertain due to the cos- 
mic variance. Moreover, we know that the Local 
Group is surrounded by the large-scale filamen- 
tary structures in the density field such as the Lo- 
cal Superclustcr, Perseus-Pisces supercluster, the 
Great Attractor, etc. They are expected to con- 
tain warm/hot gas and could serve as targets for 
its detection. It would be interesting therefore to 
explore observational signatures of this gas and 
make theoretical predictions on the sky coverage 
and column density distributions of gas in these 
structures. 

Fortunately, this problem can be addressed by 
using constrained simulations (Klypin et al. 2001a, 
hereafter Paper I), i.e., cosmological simulations 
with initial conditions constrained to reproduce 
the observed large-scale structures in the nearby 
Universe (see also Kolatt et al. 1996; Bistolas & 
Hoffman 1998; van de Weygaert & Hoffman 1999). 
The initial conditions are generated within the 
framework of a flat CDM cosmogony by means of 
constrained realizations of Gaussian fields, where 
the constraints are set by the MARK III survey of 



peculiar velocities. The setup of these realizations 
is such that that the structure on scales larger than 
a few Mpc is tightly constrained by the data, while 
smaller scales essentially constitute a random re- 
alization of the assumed primordial perturbation 
field. 

In the study presented here, we used the 
method of Klypin et al. (2001a) to set up and 
run constrained iV-body+gasdynamics simula- 
tions of the Local Superclustcr region (LSC; a 
roughly spherical region of R < 30/i _1 Mpc cen- 
tered on the Virgo cluster) in ACDM cosmol- 
ogy. The simulation followed gasdynamics of 
the baryonic component and collisionless dynam- 
ics of dark matter particles using a newly de- 
veloped eulerian adaptive mesh refinement code 
(see § 2). Gas cooling (except for adiabatic ex- 
pansion) was not included. The main results of 
our study depend on the physics of shock-heated 
(T > 10 5 K) gas in the relatively low density 
(p/(p) ~ 1 — 10) regions. The cooling time of such 
gas is i CO oi » 3.1 x 10 11 yrs K- 1 A_ 23 (1 + 5)-^, 
where we assumed fif, = 0.02ft, -2 and complete 
ionization of gas, A_ 2 3 is cooling rate in units of 
10~ 23 cm 3 erg s _1 , S is gas overdensity, and T 5 is 
gas temperature in 10 5 K. The cooling time is thus 
longer than the Hubble time in our regions of in- 
terest and effects of cooling can be neglected. The 
cooling time can be short in the cluster cores where 
gas cooling can result in a cooling flow. Such re- 
gions, however, occupy a very small fraction of the 
volume and will not affect the main conclusions 
of the paper about the probability for the WHIM 
detection via metal absorption. The cooling will 
affect our results on the X-ray fluxes from regions 
in the cluster cores; our fluxes should be viewed 
as an upper limit on the possible emission. 

The use of observational constraints dictates 
that the simulations are initialized and evolved in 
the supergalactic coordinates used to map the ac- 
tual universe. This enables us to construct sky 
maps of the desired quantities (e.g., X-ray surface 
brightness, or column densities) corresponding to 
real objects and actual locations in the real uni- 
verse. The large angular extent of the LSC and 
its favorable orientation on the sky 2 imply that a 

2 The concentration of nearby galaxies along the supergalac- 
tic plane extends over more than 100°. The plane is almost 
perpendicular to the galactic plane; its center, the Virgo 
cluster, is close to the North Galactic Pole. 
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large fraction of observational lines-of-sight in the 
northern galactic hemisphere pass through the re- 
gions of high gas density. We analyze the distri- 
bution and properties of gas in the LSC region 
focusing on its possible observational signatures 
such as absorption and emission by OVI, OVII, 
and OVIII, soft X-ray emission, and Sunyaev- 
Zcldovich effect. 

The paper is organized as follows. The numeri- 
cal code, initial conditions and simulations are de- 
scribed in § 2 and § 3, respectively. The geometry 
of density field, properties of gas in the simulation, 
sky maps of various quantities, and observational 
signatures of the gas are discussed in § 4. We 
conclude the paper with a brief discussion of the 
results and our conclusions in § 5. 

2. Numerical code 

The adaptive refinement tree (ART) code 
(Kravtsov et al. 1997; Kravtsov 1999) achieves 
high spatial resolution by adaptively refining re- 
gions of interest using an automated refinement 
algorithm. The TV-body part of the code was 
used extensively to run high-resolution dissipa- 
tionless cosmological simulations (e.g., Kravtsov 
1999, and references therein). Due to the use of 
eulerian mesh hierarchy in the ART algorithm, the 
inclusion of eulerian gasdynamics is a natural ex- 
tension. The adaptive mesh refinement methods 
for numerical hydrodynamics enjoy widespread 
popularity in the physics and engineering com- 
munities and are now also gaining popularity in 
cosmology and astrophysics (e.g., Truelove et al. 
1998; Norman & Bryan 1999; Khokhlov et al. 1999; 
Abel et al. 2000; Plewa & Miiller 2001; Refregier 
& Teyssier 2001; Rickcr 2001). The adaptive mesh 
refinement is particularly attractive for cosmology 
because interesting regions in cosmological sim- 
ulations usually occupy only a small fraction of 
the computational volume and thus can be refined 
with relatively small number of mesh cells. 

Here we will briefly describe the gasdynamics 
part of the algorithm, detailed account of which 
will be presented elsewhere (see also Kravtsov 
1999). The equations of gasdynamics and parti- 
cle motion are integrated in "supercomoving" vari- 
ables (Shandarin 1980; Martel & Shapiro 1998). 
These variables are remarkable as their use almost 
completely (completely for ideal gas with 7 = 5/3) 



eliminates explicit dependence on cosmology in 
the model equations. The equations can there- 
fore be integrated using standard gravitation and 
gasdynamics solvers with no need for additional 
coefficients and corrections. 

A simulation is usually started with a uniform 
grid covering the entire computational domain. As 
the universe evolves and expands, additional res- 
olution is required in and around dense objects. 
The ART code uses mesh refinement to increase 
spatial resolution in regions where such increase 
is needed. The refinements are recursive: the re- 
fined regions can also be refined, each subsequent 
refinement having half of the previous level's cell 
size. This creates a hierarchy of refinement meshes 
of different resolution, size, and geometry covering 
regions of interest. Because each individual cubic 
cell can be refined, the shape of the refinement 
mesh can be arbitrary and effectively match the 
geometry of the region of interest. This strategy 
is particularly well suited for simulations of a se- 
lected region within a large computational box, as 
in the constrained simulations presented below. 

During a simulation, spatial refinement is ac- 
companied by temporal refinement. Namely, each 
level of refinement, I, is integrated with its own 
time step At; = Ar /2 ( , where Ar is the global 
time step of the zeroth refinement level in the 
code's time units (the code uses time variable de- 
fined as follows: dr = a~ 2 dt/to, where t is the 
physical time, a is expansion factor, and to is an 
arbitrary time constant). This variable time step- 
ping is very important for accuracy of the results. 
As the force resolution increases, more steps are 
needed to integrate the trajectories accurately. To 
increase mass resolution in certain regions the code 
uses particles of different masses. In the simula- 
tions presented in this paper this ability was used 
to achieve high mass resolution inside a region cen- 
tered around the Virgo cluster (see below). 

The main features of the gasdynamics im- 
plementation follow closely the algorithm of 
Khokhlov (1998). A second-order Godunov-type 
solver (Colella & Glaz 1985) is used to compute 
numerical fluxes of gas variables through each cell 
interface, with "left" and "right" states estimated 
using linear piecewise linear reconstruction (van 
Leer 1979). The Riemann solver was written and 
kindly provided to us by A.M. Khokhlov and has 
good shock-capturing characteristics (Khokhlov 
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Fig. 1. — Results of the Af-body+gasdynamics simulations of the "Santa Barbara" cluster using the ART code. The 
panels show various spherically averaged profiles; counter clockwise from the upper left: dark matter density; gas 
density; gas pressure; gas entropy. The ART profiles (triangles) are compared to the corresponding mean profiles 
from the SB cluster comparison project (Frenk et al. 1999) shown by the solid line, which represent average over 
12 (13 for the pnu profile) simulations performed using different cosmological codes. The open circles show profiles 
from the adaptive mesh refinement (AMR) simulation by Greg Bryan performed as part of the SB project. The 
open triangles show profiles for a simulation with the dark matter particle of 1 x 10 9 Mq (~ 116, 000 particles within 
the cluster virial radius), while solid triangles show results of a simulation with eight times better mass resolution 
(~ 928, 000 particles within the virial radius). The profiles from the ART simulations are shown at scales larger then 
4 formal resolutions. 
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1998; Kravtsov 1999) which is desirable for cos- 
mological applications. The gasdynamics is cou- 
pled to the dynamics of dark matter through the 
common potential, which is used to compute ac- 
celerations for both the DM particles and the 
gas. The code employs the same data structures 
and similar refinement strategy as the previous 
iV-body versions of the code. However, in addi- 
tion to the DM density criteria it also allows for 
additional refinement criteria based on the local 
density of gas, shock, sharp gradient indicators. 
The refinement criteria can be combined with dif- 
ferent weights allowing for a flexible refinement 
strategy that can be tuned to the need of each 
particular simulation. 

Extensive tests of the TV-body code and com- 
parisons with other numerical TV-body codes can 
be found in Kravtsov (1999) and Knebe et al. 
(2000). To compare performance of the N- 
body+gasdynamics code to other existing cos- 
mological gasdynamics codes in a realistic cosmo- 
logical setting, we have simulated the cluster used 
in the "Santa Barbara (SB) cluster comparison 
project" (Frenk ct al. 1999). This project com- 
pared simulations of the same Coma-size cluster 
performed using 12 different gasdynamics codes. 
Figure 1 shows spherically averaged profiles of 
dark matter density, gas density, gas pressure, and 
gas entropy for the ART simulations with two dif- 
ferent mass resolutions for the dark matter. The 
ART results are compared to the average pro- 
files from the SB project and to the profile of the 
Adaptive Mesh Refinement (AMR) simulation 3 by 
G.L. Bryan performed as part of the SB project. 
The refinements in the ART runs were based on 
the local gas density, with the threshold density 
for refinement level / equal to 4 x (m gas ) Axq 31 , 
where (m gas ) and x are average gas mass per 
and size of zeroth-level grid cell, respectively. The 
simulations used 128 3 uniform grid and six levels 
of refinement, which corresponds to the smallest 
cell size of « 8 kpc in the 64 Mpc computational 
box which is similar to the formal resolution of 
the AMR simulation. 

The figure shows very good general agreement 
between the ART and Bryan's AMR profiles, es- 

3 Note that this is entirely independent AMR code that uses 
different solvers for both gas and dark matter and different 
refinement algorithm. 



pecially for the lower-resolution ART run which 
has the same mass resolution as the AMR run. 
The higher mass resolution run produces a some- 
what steeper dark matter density profile in the 
innermost regions. This is in agreement with re- 
cent convergence studies for dark matter density 
profiles (Ghigna et al. 2000; Klypin ct al. 2001b). 
It is interesting to note that the steepening of the 
DM profile affects gas profiles, resulting in a some- 
what higher central gas density (and smaller core 
radius) and higher central temperature and pres- 
sure. The entropy profiles of the two eulerian mesh 
refinement simulations agree very well and show 
presence of a well resolved core. There is a clear 
and alarming mismatch of these profiles and the 
average SB profile which is dominated by the SPH 
simulations at small scales. The origin of the mis- 
match is unclear, although the fact that the two 
different eulerian codes agree well probably indi- 
cates that it is due to less accurate treatment of 
shocks in the SPH codes. We are currently work- 
ing on more detailed comparisons of entropy evo- 
lution in ART and SPH codes. 

3. Constrained simulations 

The goal of this paper is to perform numeri- 
cal simulations that match the observed local uni- 
verse as well as possible. Namely, we are inter- 
ested in reproducing the observed structures: the 
Virgo cluster, the Local Supercluster (LSC) and 
the Local Group (LG), in the approximately cor- 
rect locations and embedded within the observed 
large-scale configuration dominated by the Great 
Attractor and the Perseus-Pisces supercluster. 

The setup details of the constrained initial con- 
ditions used for our simulations arc presented in 
(Klypin et al. 2001a). Initial conditions arc con- 
structed using constrained realizations based on 
the MARK III catalog of radial peculiar velocities 
and assuming the ACDM model with the follow- 
ing parameters: flo = 0.3, Cl\ = 0.7, as = 0.9, 
^baryon = 0.03, h = 0.7. Constrained realiza- 
tions of Gaussian fields are generated by sampling 
the random residual from the mean (Wiener fil- 
ter) field given the data and the assumed model 
(Hoffman & Ribak 1991). The nature of the data 
and the ACDM model implies that scales larger 
than a few Mpc are tightly constrained by the 
data while the small scale structure constitutes 
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a random realization of the assumed model (see 
Klypin et al. (2001a) for a quantitative analysis). 
The random (unconstrained) waves dominate on 
small scales which leads to random displacements 
of some large-scale features in the density distri- 
bution by a few Mpc from realization to realiza- 
tion. The realization used for our simulations was 
selected from a set of low-resolution dissipation- 
less runs to have the position of the Virgo cluster 
to be as close to its observed location as possible. 
The realizations are constructed assuming that the 
linear theory is valid on all scales at the present 
epoch, and is then evolved backwards in time to 
the initial epoch of the simulation by the linear 
growth factor. 

A region within a radius of 30/i _1 Mpc around 
the Virgo cluster inside the 160/i _1 Mpc computa- 
tional box was selected for re-simulation at the 
highest mass and force resolution. The initial 
positions and velocities for 1024 3 particles were 
generated in the entire 160/i _1 Mpc box (roughly 
centered on the Virgo cluster). Particles outside 
the designated high-resolution region were then 
merged 4 into particles of larger mass and this pro- 
cess was repeated for merged particles to con- 
struct a nested hierarchy of particles of four differ- 
ent masses surrounding the high-resolution region. 
The simulation used four DM particle species with 
DM particle mass of 3.16 x 1O 8 /i _1 M in the high- 
est resolution region. The initial conditions used 
for our simulations are identical to the simulation 
analyzed in Klypin ct al. (2001a). 

The simulations presented here were run using 
128 3 zeroth-level grid. The initial gas density and 
velocity for each grid cell was computed using DM 
particle positions and velocities and cloud-in-cell 
(CIC) interpolation. The subsequent mesh refine- 
ment was done using criteria based on the local 
DM and gas densities. Namely, a cell of level 
L is marked for refinement if mass of DM par- 
ticles (evaluated using the CIC-smoothed density 
field) and/or mass of gas in the cell exceeds spec- 
ified threshold. The final refinement map is then 
created by smoothing the original map using an 
algorithm similar to that described in Khokhlov 
(1998). The smoothing generally enlarges regions 

4 The larger mass (merged) particle is assigned a velocity 
and displacement equal to the average velocity and dis- 
placement of the smaller-mass particles. 



marked for refinement and reduces noisiness of 
the original map. The threshold can be chosen 
to be different at different levels of refinement. 
The density-based mesh refinement is appropriate 
as we are interested in resolving collapsing high- 
density regions. However, the specific choice of 
threshold masses is arbitrary. 

To ensure that the results of our study do 
not depend on the choice of thresholds, we have 
run two simulations with very different numeri- 
cal criteria for refinement. In the first simulation 
(hereafter Rl), the mass thresholds were the same 
for all levels and were set to the average mass 
of DM and gas per zeroth level cell ((todm) ~ 
1.4 x IO^^Mq and (m g ) w 2.3 x IO^/j^M©, 
respectively). There was no limit on the num- 
ber of refinement levels; a total of 8 levels were 
introduced in the highest density regions, corre- 
sponding to the smallest cell size of ~ 4.9/i _1 kpc. 
The second simulation (hereafter R2) was started 
from the same initial conditions and used the same 
threshold masses for refinement on the zeroth and 
the first levels as the simulation Rl. However, run 
R2 was limited to four levels of refinement (cor- 
responding to the smallest cell size of 78/i -1 kpc 
and mass thresholds for levels 2 and 3 were set 
to (m)/4 and (m)/16, respectively. The simula- 
tion R2, therefore, has lower peak resolution than 
the Rl. The fraction of the volume resolved at 
the peak resolution is, however, much greater in 
R2 than in Rl. In particular, all regions of local 
overdensity ~ 20 — 30 were refined to the peak res- 
olution of 78/i -1 kpc in R2. In the simulation Rl 
this resolution was reached only in regions with 
local overdensity higher than ~ 500. For all the 
results presented in this paper there is no signifi- 
cant difference between simulations Rl and R2. In 
the remainder of the paper we will use simulation 
R2 due to its higher resolution in the intermediate 
density regions. 

The time step on the zeroth level grid for the 
simulations, A To, was set at the beginning of each 
global step to satisfy the Courant condition (with 
the Courant-Friedrichs-Lewy factor of cfl — 0.4) 
for cells on all levels and with the condition that 
the corresponding step in expansion factor would 
not exceed Aa < 0.003 (the latter condition is 
usually imposed only during the early stages of 
evolution). The time steps for particles and cells 
on other levels were scaled appropriately: At; = 
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Atq2~ 1 . The simulations were started at z = 30 
and a total of 364 global time steps were performed 
in both simulations. 

4. Results 

4.1. Spatial distribution of gas and dark 
matter 

At the current epoch most of the mass (« 
7.5 x 10 14 /i- 1 M Q ) of the LSC is located in a fil- 
ament roughly centered on the Virgo cluster and 
extending over ~ 40/i _1 Mpc. The simulated Local 
Group (LG) is located in an adjacent smaller fila- 
ment, which is not a part of the main body of the 
LSC, and has a peculiar velocity of w 250 km s _1 
toward the Virgo cluster. The peculiar velocity 
flow in the vicinity of the LG in the simulation 
is relatively "cold": the peculiar line-of-sight ve- 
locity dispersion within 7/i~ 1 Mpc of the LG is 
< 60 km s _1 , comparable to the observed velocity 
dispersion of nearby galaxies. In the main filament 
of the LSC, the peculiar velocities are higher with 
typical values of ~ 200 — 300 km s _1 ; the peculiar 
velocities are the highest, > 400 — 500 km s _1 in 
the immediate vicinity of the Virgo cluster. The 
overall matter distribution and peculiar velocity 
field in the the simulated LSC region is discussed 
in detail in Paper I and we refer reader to this pa- 
per for further information. Here we will briefly 
discuss the distribution of gas and dark matter. 

Figure 2 shows a 15/i _1 Mpc slice through the 
gas and DM distribution in the simulation R2. 
The slice is centered on the simulated Virgo clus- 
ter. The left panels show the (continuous) density 
field of gas (upper left panel) and distribution of 
DM particles (lower left panel). Comparison of the 
gas and DM distributions shows that at these large 
scales gas traces the gravitationally dominant DM 
distribution very well: similar filamentary struc- 
tures and halos can be identified in both panels. 
As discussed above, the dominant structures are 
the "Virgo" cluster in the center and the hori- 
zontal filament around it. The right panels show 
the projected emission-weighted temperature field 
of the gas (upper panel) and projected density of 
the Warm-Hot (solid contours) and diffuse (dotted 
contour) phases of the gas. The panels show that 
warm (~ 10 5 — 10 6 K) gas is located in filaments, 
while hotter gas (> 5 x 10 6 K) is concentrated 
around virialized halos of groups and clusters. It 



is the Warm-Hot gas that will be the main focus 
of our study. 

Below we discuss properties and possible obser- 
vational signatures of the intergalactic gas in the 
constrained simulations of the LSC region. The 
analysis presented below excludes the inner struc- 
ture of the Virgo cluster. The cooling and feed- 
back processes neglected in our simulations would 
affects the structure of the cluster. We will present 
a constrained simulation of the Virgo cluster that 
includes heating and cooling processes in a forth- 
coming paper. 

4.2. Properties of gas 

Dave et al. (2001) present a detailed discussion 
of gas properties in the currently favored struc- 
ture formation models. They classify the gas in 
four phases: diffuse (d < 1000, T < 10 5 K), con- 
densed (S > 1000, T < 10 5 K), hot (T > 10 7 K), 
and warm- hot (10 5 < T < 10 7 K), where 5 and 
T are gas overdensity and temperature, respec- 
tively. The simulations presented here do not in- 
clude cooling and thus cannot predict the frac- 
tion of gas in condensed phase, while fraction of 
gas in the hot phase is sensitive to the presence 
or absence of a small number of galaxy clusters 
in the simulation volume. However, the hot gas 
fraction is expected to be small (~ 5%) and the 
gas not in diffuse or WH phase should correspond 
approximately to the condensed fraction in the 
simulations with cooling. It is therefore inter- 
esting to compare the fraction of WH and dif- 
fuse gas to those in previous simulations. Dave 
ct al. (2001) find that - 40% and - 30% of the 
gas at present epoch is in the diffuse and Warm- 
Hot (WH) phases, respectively. The fraction of 
WH gas appears to be rather robust: approxi- 
mately the same fraction of WH gas is measured 
in simulations with and without cooling and/or 
supernova feedback. The corresponding fractions 
of diffuse and WH gas in the high-resolution re- 
gion of our simulations are Sldiff/^ ~ 0.41 and 
flwH/^b ~ 0.34, respectively. Both fractions are 
in good agreement with the fractions found by 
Dave et al. (2001). 

Figure 3 shows the X-ray volume emissivity in 
the [0.5 — 2] keV energy range and gas temperature 
as a function of gas density in the simulation Rl. 
The points in the figure represent individual mesh 
cells randomly selected from the high-resolution 
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Fig. 2. — A slice through the gas and DM distribution in the simulation R2. The slice is centered on the simulated 
Virgo cluster and is 15h~ Mpc thick. Upper left panel: the contour plot of the projected density field, where the 
thick solid contours correspond to the mean gas density, while subsequent contours correspond to the densities of 
4, 8, 16, etc. times the mean density. Lower left panel: the distribution of DM particles in the slice; randomly 
selected subset of 10% of all the particles in the high-resolution region is shown. Upper right panel: the contour 
plot of the emission-weighted projected temperature field in the slice. The thick contour corresponds to temperature 
of 3 x 10 5 K, while other contours correspond to temperatures of 10 6 K, 4 x 10 6 K, 6 x 10 6 K, 8 x 10 6 K, and 
10 7 K. Lower right panel: the projected density fields of the Warm-Hot (solid contrours) and diffuse phases of the 
gas (dotted contour). The thick solid contour corresponds to the density of WHIM equal to the mean density of the 
gas, while thin solid contours correspond to the WHIM densities of 8, 16, 32, etc. times the mean gas density; the 
dotted contour corresponds to the density of diffuse gas equal to 0.1 of the mean gas density. 
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Fig. 3. — X-ray (0.5-2 keV) emissivity (top panel) 
and gas temperature (bottom panel) as a function of 
gas density. The gas density is given in units of the 
mean density. The X-ray emissivity was computed 
using the Raymond & Smith (1977) code assuming 
uniform metallicity of 0.3 solar. Points correspond to 
randomly selected individual mesh cells in the high- 
resolution simulation region. Note that the cells plot- 
ted correspond to relatively high-density regions and 
represent less than half of the volume; the median val- 
ues of gas density, temperature, and X-ray emissivity 
by volume (i.e., half of the volume is at values below 
the median) are 0.29, 5000 K (the minimum assumed 
gas temperature), and 10 -42 ergs s _1 cm -3 . 

region of the simulation (see § 3). The X-ray emis- 
sivity was calculated using the Raymond & Smith 
(1977) code assuming complete ionization and uni- 
form gas metallicity of 0.3 solar. Figure 3 shows 
that temperature is correlated with gas density, 
although there is much scatter (especially at low 
gas densities). As noted by Dave et al. (2001), the 
average can be roughly described as p oc T. The 
large scatter in temperature translates into signif- 
icant scatter in X-ray emissivities. The relatively 



sharp upper edge in the ex — Pg scatter plot is due 
to the fact that at high tmperatures (correspond- 
ing to the highest emissivities), the emissivity is 
a strong function of gas density (oc p g ) and is a 
relatively weak function of temperature (oc T 1 / 2 ). 
As can be seen from the lower panel, the maxi- 
mum temperatures of the gas change by less than 
a factor of ten over four orders of magnitude in 
density. The upper envelope in the upper panel is 
therefore dominated by the density dependence: 
e x oc p 2 g and scatter in temperature intoduces only 
a small scatter at the maximum emissivity at a 
given p g . 

Note that although a large fraction of WH gas 
is at relatively low overdensities (p g /(p g ) < 100), 
the X-ray emission is dominated by gas at higher 
overdensities. Note also that at overdensities of 
Pg/(p g ) ~ 10 most of the gas has temperatures 
in the range T ~ 10 5 - 10 6 K. The correlation of 
emissivity and temperature with density presented 
in Figure 3 are in good agreement with results 
of previous studies (Dave et al. 2001; Croft et al. 
2001, see their Fig. 6 and Fig. 2, respectively). 

4.3. Sky maps 

The use of observational constraints implies 
that the computational box axis in our simula- 
tions are mapped with the same supergalactic 
coordinates as the observed universe. Because ob- 
servational constraints small scales were not too 
tight, there was a room for fine-tuning the position 
of the origin of coordinates (i.e., position of the 
Local Group in the simulations). We choose the 
origin to be the location of a pair of galaxy size 
halos in the high-resolution dissipationless sim- 
ulation described in Paper I (see Fig. 6 in that 
paper). For this choice of origin, the simulated 
Virgo cluster has supergalactic coordinates of 
(SGX,SGY,SGZ) = -1.2, 13.7, OJ/i^Mpc and 
is located 13.74/i _1 Mpc from the Local Group. 
Knowledge of the supergalactic coordinates allows 
us to construct sky maps in real celestial coordi- 
nates. 

In the following we present a series of sky maps 
showing projections of various quantities in a re- 
gion of the sky around the north galactic pole. 
The maps are shown as polar views in the galactic 
coordinates, with the North Galactic Pole (NGP; 
b = 90°) in the center. We plot only a region 
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so 



90 




of b > 30°. We do not consider other regions of 
sky because 1) most of the nearby structures in 
the high-resolution region (w 30ft. _1 Mpc around 
the Virgo cluster) of the simulations are concen- 



('-^gHrated around the NGP (the Virgo cluster is lo- 
3.0 cated within 15° of the NGP) 2) the sky at low 
2 j galactic latitudes is more difficult to observe due 
to galactic obscuration. All maps include only 
I cells in the high-resolution region of the simula- 
Bo,5 tion; the resolution outside this region was con- 
H_ g siderably worse (cell size of 1.25ft, _1 Mpc) and our 
H predictions would be considerably less accurate. 

1 - 1 All quantities represent averages over a number 
B-3.0 of random lines of sight passing through each map 
pixel. The integration along lines of sight was done 
from 2/i _1 Mpc to 25/i _1 Mpc. We thus exclude the 
regions in the immediate vicinity of the LG which 
are not modelled accurately in our simulations. 

Figure 4 shows the sky maps of the gas over- 
density and X-ray emission weighted temperature. 
Note that although the gas distribution is very fil- 
amentary in appearance, the overall distribution 

lcg(T x ) 

is flattened along a plane which is roughly aligned 
7.0with the observed supergalactic plane. The large 
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Fig. 4. — Sky maps of density (top panel) and X- 
ray emission weighted temperature (bottom panel) of 
the gas in the simulated Local Supercluster. The pan- 
els show polar views around the North Galactic Pole 
(b > 30°). Radial lines and circles correspond to spac- 
ings of 30° and 15° in galactic / and b, respectively; 
I = 0° is at the bottom and I = 90° is on the left; 
b — 90° is in the center. The density represents the 
average maximum overdensity along a line of sight in 
a given map "pixel"; the temperature is in degrees 
Kelvin and is avaluated as an emission-weighted av- 
erage along lines of sight in a "pixel" spanning from 
2/i _1 Mpc to 25/i~ 1 Mpc from the origin. The corre- 
sponding color scales are shown by the bar to the right 
of each panel. The Virgo Cluster is close to the cen- 
ter of the plots; the Ursa Major cluster is close to the 
top of the plots. The Local Void is in the bottom left 
quadrant. 



g e v °id in the lower left corner is the counterpart 
of the Local Void in the observed distribution of 
6 -%iearby galaxies. The temperature of gas in fila- 
6 [ments ranges from one to a few million degrees K 
and reaches tens of million K within and around 
5 ' 6 the virialized regions of groups and clusters. Note 
5.3that these regions are surrounded by the narrow 
regions of enhanced temperature corresponding to 
accretion shocks. The accretion shocks are also 
visible around massive filaments. The two hottest 
regions in the center and at the top of the plot 
are the environments of the Virgo and Ursa Ma- 
jor clusters, respectively. 

Figure 5 shows sky maps of surface brightness 
of X-ray emission in [0.5 — 2] keV energy band and 
column densities of oxygen ions: OVI, OVII, and 
OVIII. The hot, T > 10 6 K gas in filaments and 
virialized regions (see Figures 3 and 4) should 
emit soft X-ray radiation. As can be seen from 
the Figure 5, the X-ray surface brightness map 
has a more patchy appearance than both the den- 
sity and temperature maps. This is because X-ray 
emission is roughly proportional to p 2 T 5 and is 
thus heavily weighted towards the highest density 
regions. 

The emission- weighted average X-ray spectrum 
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from the considered sky region is shown in Fig- 
ure 6. The spectrum was computed using the 
Raymond & Smith (1977) plasma code and as- 
suming uniform gas metallicity of 0.3 solar. The 
model spectrum is compared to the broad-band 
measurements of soft X-ray extragalactic back- 
ground from the RASS survey with the contribu- 
tion of AGNs, stars, and the Local Hot Bubble 
subtracted (see Kuntz et al. 2001, for details). In 
both models and observations emission from the 
3° region around the Virgo cluster was excluded. 
The figure shows that energies < 0.5 keV the pre- 
dicted emission of the intergalactic gas is signifi- 
cantly lower than the measured background. As 
argued by Kuntz et al. (2001), this emission is 
likely due to the Milky Way's own hot halo. How- 
ever, at energies ~ 1 keV the X-ray emission of 
the intergalactic gas, although still relatively small 
(~ 10% of observed background at these ener- 
gies; the measured total background at 1 keV is 
~ 10 keV s _1 cm _2 sr _1 keV _1 ), may be important 
for understanding the total background at these 
energies. 

In addition to the X-ray brightness map, Fig- 
ure 5 shows the sky maps of column densities of 
three ionic species of oxygen: OVI, OVII, and 
OVIII. The concentrations of these species were 
calculated using the current version (94.00) of the 
CLOUDY code (Ferland 2000) assuming that the 
gas is embedded in the observed XRB spectrum 
and has uniform metallicity of 0.3 solar. We as- 
sume that power-law XRB spectrum. P(E) = 
AE~ r , with A = 9.5 keV s^cm^sr^keV- 1 
(Kuntz et al. 2001) and T = 1.46 (Chen et al. 
1997); this is equivalent to the flux J„ = 1.2 x 
10~ 26 (zV^x)~ ' 46 er S s s _1 cm~ 2 sr _1 Hz _1 , where 
hvx = 40 keV. We also assume that the spec- 
trum has an exponentials cutoff at energies above 
hvx and at very small energies (<C 1 keV). 

Comparison of the column density maps in 
Figure 5 to the projected density map in Fig- 
ure 4 shows that high column density regions are 
roughly correlated with regions of high gas density. 
Note, in particular, that high-column density re- 
gions are in general located near the supergalactic 
plane. Nevertheless, different ionic species trace 
regions of somewhat different densities and tem- 
peratures (e.g. Hellsten et al. 1998; Nahar 1999) 
and their maps, therefore, differ in details. For 
instance, the column densities of OVI and OVII 



a 



io- 3 



- 


□ 


□ 


iii 


b > 30° 


- 








□ 






E 
- 






□ 




: 
- 




RASS 












central 


3" of Vir| 


ro excluded 



















0.1 1 
E [keV] 



Fig. 6. — The X-ray spectrum of gas in the LSC re- 
gion. The spectrum, shown by the green solid line, 
represents the emission-weighted average over the sky 
region around the NGP (6 > 30°) shown in Figure 5. 
Open squares show the broad-band measurements of 
the intensity of the X-ray background by the RASS 
with known contributions of AGNs, stars, and the Lo- 
cal Hot Bubble subtracted (Kuntz et al. 2001). The 
3° region around the center of the Virgo cluster was 
excluded in both model and observations. 

are depressed in the vicinity of the Virgo cluster 
(I ~ 300°, 6- 85 -90°) due to the high gas tem- 
peratures around the Virgo cluster. The decrease 
in column density of OVIII in this region is rela- 
tively small. 

Our primary interest is evaluating the prospects 
of detection of the gas in the LSC region. The oxy- 
gen is the most abundant of the metals that can 
absorb or emit in UV and X-rays. The absorption 
or emission by lines of the ionic species of oxy- 
gen is therefore the most promising probe of the 
LSC gas. The strongest lines of OVII and OVIII 
are in soft X-rays at 0.5740 keV and 0.6536 keV, 
respectively. The OVI can be detected using reso- 
nance line doublet in UV at AA1031.92, 1037.62 A. 
The current instruments in X-ray (Chandra and 
XMM-Newton) and UV (FUSE) are sufficiently 
sensitive to detect at least some of the intergalac- 
tic gas. The exact detection limit depends on 
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Fig. 5. — Sky maps of X-ray surface brightness in [0.5 — 2] keV (top left) and column densities of oxygen ions: OVI 
(top right), OVII (bottom left), and OVIII (bottom right). The X-ray brightness is given in keVs _1 cm _2 sr _1 keV _1 , 
while column densities are in cm -2 . The layout is similar to that of Figure 4. The brightness and column densities 
were evaluated by averaging along random lines of sight within each map "pixel" in the distance range [2, 25]ft- 1 Mpc. 
See text for details. 
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the brightness of the background source (in the 
case of absorption) and specifics of the instru- 
ment and observation (time of integration, spec- 
tral resolution, etc.). For an X-ray bright quasar 
and a reasonable integration time (100 ks), the 
Fang & Canizares (2000) estimate that Chandra, 
XMM, and the planned Constellation-X should de- 
tect OVIII (sensitivity to OVII is similar) at col- 
umn densities higher than 3.5 x 10 16 , 7.6 x 10 15 , 
and 1.3 x 10 15 cm -2 , respectively. The OVI was 
robustly detected in absorption with FUSE for col- 
umn densities > 10 14 cm~ 2 (Tripp et al. 2000). 

The column density maps in Figure 5 show that 
for the assumed metallicity the column densities of 
OVI should be detectable with FUSE over a large 
fraction of the sky, especially near the supergalac- 
tic plane. On the other hand, the column densi- 
ties of OVII and OVIII detectable with Chandra 
and XMM-Newton cover a much smaller fraction 
of the sky. This fraction is much higher for the 
sensitivity of the Constellation-X. 

The cumulative fraction of the sky (estimated 
over the region of the sky shown in Figure 5, i.e. at 
b > 30°) at column densities above a given value is 
shown in Figure 7. This fraction can be thought of 
as a probability to have a detectable column den- 
sity in a random line of sight. The figure shows 
predictions under three different models for gas 
metallicity: uniform metallicity of 0.05 and 0.3 so- 
lar, and density-dependent law, 

Z/Z Q = 0.025(p 5 /<p fl )) - 35 , (1) 

which approximates results of numerical simula- 
tions of Cen & Ostriker (1999a). For uniform 
metallicity of Z/Zq = 0.3, the probability is ~ 0.5 
for the detection of OVI absorption with FUSE. 
The OVIII absorption (results for OVII are quanti- 
tatively similar) should be detected with the prob- 
ability of « 0.5 with Constellation-X and < 0.05 
with the XMM-JVewton and Chandra. For lower 
metallicities of the gas, the probabilities are con- 
siderably lower. Note that sky fraction distri- 
butions for the density-dependent metallicity are 
similar to the uniform metallicity of Z/Zq = 0.05. 
This indicates that the signal is dominated by gas 
at overdensities of ~ 5 — 10 (see eq. 1) and is sen- 
sitive to the average metallicity of this gas. Note 
that the fractions are calculated over the entire sky 
region around the NGP (b > 30°), which includes 
sizeable voids in the gas density; the fractions will 




Fig. 7. — The fraction of the sky around the NGP 
(b > 30°) at column densities greater than a given 
value N for OVI and OVIII for different metallicity 
models. The results for OVII are quantitatively simi- 
lar to those for OVIII. The solid and dashed lines show 
results for the uniform metallicity of gas of 0.3 and 0.05 
solar, respectively; the dot-dashed shows results for 
the density-dependent metallicity of the gas: Z/Zq = 
0.025(5° 35 , where 5 g is gas overdensity, which approx- 
imates results of numerical simulations of Cen & Os- 
triker (1999a). The sensitivity of current and future 
instruments is also shown (sensitivity of the XMM- 
Newton is intermediate to those of Constellation-X 
and Chandra). Note that for the metallicity of the 
IGM of 0.3 solar, OVI and OVII/OVIII have column 
densities detectable with the FUSE and Constellation- 
X over a large fraction of this sky region. 

be higher if one focuses on the sky regions near 
the supergalactic plane. 

5. Discussion and conclusions 

We presented results of high-resolution gas 
dynamics+V-body cosmological simulations of 
the Local Supercluster region constrained to re- 
produce the observed nearby large-scale density 
field. The simulations focused on the region of 
w 30ft, Mpc around the Virgo cluster and fol- 
lowed the dynamics of dissipationless dark matter 
and adiabatic dynamics of gas using the new ver- 
sion of the Adaptive Refinement Tree code that 
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includes gasdynamics. The detailed discussion 
of constrained simulations, setup of initial condi- 
tions, and properties of the density and velocity 
fields in the simulated Local Supercluster are pre- 
sented in a separate paper (Klypin et al. 2001a). 
In this paper we use the simulations to study the 
properties and possible observational signatures 
of intergalactic medium in the LSC region. We 
find that within 30/i _1 Mpc of the Virgo cluster 
~ 30% of baryons is in the warm/hot phase at 
T ~ 10 5 - 10 7 K and « 40% are in the diffuse 
phase (T < 10 5 K, 5 < 1000). The latter phase 
occupies most (> 90%) of the volume. 

We presented simulated sky maps of 0.5 — 2 keV 
X-ray surface brightness and column densities of 
three ionic species of oxygen (OVI, OVII, and 
OVIII). The maps are constructed in real galactic 
coordinates for the sky region around the North 
Galactic Pole (6 > 30°). The map of X-ray 
brightness has a patchy appearance with most of 
the emission coming from the high-density envi- 
ronments of groups and clusters. We compared 
emission-weighted average X-ray spectrum from 
the warm/hot gas in the simulation (Figure 6) 
to the measurements of soft X-ray background 
(Kuntz et al. 2001). The predicted flux at en- 
ergies < 0.5 keV is 10 to 100 times lower than 
the observed flux. The LSC gas therefore cannot 
explain this soft X-ray emission which is likely to 
be associated with the hot halo of the Milky Way. 
At energies around 1 keV, the predicted flux con- 
stitutes w 5 — 10% of the total XRB flux and is 
therefore an important component of the XRB at 
these energies. At higher energies, the contribu- 
tion of the LSC gas to the XRB is insignificant. 
It is not clear whether the X-ray emission of the 
intergalactic gas can be reliably detected in the 
near future. Such detection is difficult because one 
needs a survey that covers a large sky area and is 
sufficiently sensitive to detect the signal only ~ 1% 
of the XRB. The diffuse LSC emission should be 
correlated with the supergalactic plane, but so is 
the X-ray emission from nearby AGNs (Shaver & 
Pierre 1989). 

Recently, Boughn (1999) analyzed the HEAO 
1 full sky 2 — 10 keV X-ray map and found ev- 
idence for the "diffuse" (unresolved) X-ray emis- 
sion associated with the supergalactic plane. The 
quoted maximum surface brightness of the diffuse 
emission constitutes about 1% of the 2 — 10 keV 



XRB and is Ix ~ 5 x 10~ 10 ergs s _1 cm _2 sr _1 
(see his Table 1). This corresponds to ~ 0.04 keV 
s _1 cm _2 sr _1 keV _1 assuming that flux is constant 
over 2 — 10 keV. This flux is consistent with the 
predicted flux from the WH gas in our simulation 
(see Fig. 6): 10~ 2 - lO^keV s^cm^sr^kcV- 1 
at energies 2 — 5 keV (i.e., w 0.1 — 1% of the total 
XRB at these energies; see Kuntz et al. 2001); the 
flux decreases steeply at higher energies. 

As can be seen from the X-ray brightness map 
in Figure 5, most of the X-ray emission is indeed 
concentrated towards the supergalactic plane. 
However, the emission is patchy and is far from 
being uniform. In contrast, Boughn (1999) mod- 
elled the gas distribution using a simple "pillbox" 
model for the distribution of gas in the LSC re- 
gion: the uniform gas distribution within a disk 
of radius i?sc and thickness H$c ~ 0.25i?sc- For 
the detected diffuse X-ray flux this model implies 
gas density of 2.5 x 10~ 6 cm~ 3 for temperatures of 
around 10 keV. This temperature is much higher 
than the typical temperatures of the LSC gas 
in our simulation. The difference is due to the 
fact that the gas distribution in the simulated 
LSC is not described by the pillbox model. As 
can be seen from Figures 4 and 5, the distribu- 
tion of matter in the LSC region is filamentary 
rather than disk-like and the X-ray emission is far 
from being uniform. The bulk of the X-ray emis- 
sion thus comes from the relatively high-density, 
n e ~ 10~ 5 — 10~ 3 cm~ 3 , regions within and around 
groups and clusters. 

The hot and dense regions of the LSC should 
also distort the cosmic microwave background ra- 
diation via the inverse Compton or Doppler scat- 
tering, the thermal and kinetic Sunyaev-Zel'dovich 
effect, respectively (SZ; Sunyaev & Zeldovich 
1980, and references therein). The temperature 
fluctuations of the CMB due to the non-relativistic 
thermal SZ effect can be written as: 

AT (e x + l \ ... 

where x = hv/kT CMB , y = 1.117 x 10~ 34 / n e T e dl 
(all quantities are in cgs units and the integral is 
along linc-of-sight). In the Raylcigh- Jeans regime 
(hv <C fcTcMB): AT/T w 2y; the deviation of 
CMB temperature along a given direction is thus 
proportional to the gas pressure integrated along 
this direction. 
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Although for temperatures and densities typical 
for the gas in the LSC the SZ effect is not expected 
to be very strong, due to the large angular extent 
of the LSC, it should contribute to the anisotropy 
on the on the angular scales of <~ 1 — 10°. The 
thermal SZ map has patchy appearance similar 
to that of the X-ray brightness in Figure 5 with 
AT/T - 5 x 1(T 6 - 1(T 5 within groups and clus- 
ters, <~ 10~ 7 on their outskirts, and <~ 1CP 8 in the 
strongest filaments. The kinetic SZ effect has sim- 
ilar magnitude. These fluctuations are small and 
are below current sensitivity limits of the SZ ob- 
servations; their contribution to the large angular 
scale anisotropies measured by COBE satellite is 
insignificant. 

Finally, we presented (Figure 5) sky maps of 
column densities of the three ionic species of oxy- 
gen: OVI, OVII, and OV11L The column densities 
were calculated using the densities and tempera- 
tures of gas in the high-resolution region of the 
simulation and assuming the uniform metallicity 
of 0.3 solar and the observed X-ray background. 
Although, the OVI is the least abundant (its abun- 
dance is more than an order of magnitude smaller 
than abundances of OVII and OVIII) of the three 
ions, absorption or emission in its resonant doublet 
at AA1031.92, 1037.62 A offers currently the best 
prospect for detection of the warm/hot intergalac- 
tic gas in the LSC region. For the assumed metal- 
licity of the IGM, our simulations predict that the 
probability to have a column density of OVI de- 
tectable with FUSE along a random line of sight 
in the northern galactic hemisphere is > 0.5; the 
probability is higher in the regions of the sky close 
to the supergalactic plane. 

Our predictions would be affected if metallicity 
in filaments is significantly smaller than 0.3 solar 
(the metallicities observed in the virialized regions 
of clusters). However, there are good reasons to 
believe that the filaments should be enriched to 
the levels comparable to the intracluster medium. 
The observations of metals in the Lyman alpha 
forest at high redshifts indicates that regions with 
overdensities of ~ 10 (corresponding to filaments) 
have metallicities in the range 10~ 3 — 10~ 2 solar 
(see, e.g., Cowie & Songaila 1998; Schaye et al. 
2000). The simulations predict a steep depen- 
dence of metallicity on the density in this regime 
with metallicity quickly reaching values of ~ 0.1 
solar in the vicinity of the virialized regions (e.g., 



Gnedin 1998; Cen & Ostriker 1999a; Aguirre et al. 
2001). As the Universe evolves, the metallicity 
in the low-density regions is expected to increase 
steadily (Cen & Ostriker 1999a). At z = 0, numer- 
ical simulations of Cen & Ostriker (1999a) predict 
the average metallicity of ~ 0.05 — 0.1 solar (with 
a significant scatter around the average) for over- 
densities of ~ 10; their results can be approxi- 
mated by the density-dependent metallicity given 
by equation 1. Figure 7 shows that our results 
are sensitive not to the particular density depen- 
dence of metallicity but to the average metallic- 
ity of gas at overdensities w 1 — 10. Detections 
(or lack thereof) of oxygen absorption within the 
Local Supercluster should therefore provide use- 
ful constraints on the z = metallicity of IGM in 
filaments. 

Although the structures in the simulation re- 
produce the large-scale density field rather well, 
we expect certain differences with the real ob- 
served structures. For example, the simulated 
Virgo cluster is displaced by several degrees from 
the location of the real Virgo. We cannot thus reli- 
ably predict the column density in any specific di- 
rection. However, viewed statistically, our results 
suggest that for a sample of background sources 
distributed over the northern sky, about half of 
the lines of sight should pass through the column 
density > 10 14 cm~ 2 . Unfortunately, the prob- 
ability decreases steeply for higher column densi- 
ties (w 10% for N vi > 4x 10 14 cm" 2 . We predict 
that OVI column density is very low in the area 
of w 10° around the Virgo cluster. This predic- 
tion is only for the IGM; the halos of individual 
galaxies in or around Virgo may cause OVI ab- 
sorption. The OVI absorption of the intergalactic 
gas should be easy to separate from the local ab- 
sorption in the Milky Way: most of the absorbing 
gas should be redshifted by - 3 ± 2 x 10~ 3 . The 
typical absorption should also be fairly broad as 
the typical line of sight peculiar velocities of gas in 
the high density regions are ~ 100 — 400 km s _1 
(see Figs. 5 and 6 in Paper I). 

The absorption lines of OVII and OVIII in soft 
X-rays will be difficult to detect with current in- 
struments. The predicted probability for absorp- 
tion detectable with Chandra or XMM-Newton is 
< 5% for a random line of sight in the north- 
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ern galactic sky. For the planned Constellation-X 5 
mission the probability is an order of magnitude 
higher. The simulations predict that regions of the 
highest column densities of OVII and OVIII cor- 
respond to the high-density regions in and around 
galaxy groups (see Fig. 5). 

We would like to point out that our predictions 
are made under the assumptions negligible cool- 
ing and absence of significant heating of the in- 
tergalactic gas. We think that effects of cooling 
should not have a significant effect on our predic- 
tions for the oxygen ion column densities. The 
cooling time is sufficiently short only in the high 
density regions which would occupy a very small 
fraction of the sky and thus would not affect the 
cumulative distribution of column densities. Com- 
parison of adiabatic simulation to the simulations 
that include cooling and stellar feedback presented 
by (Dave et al. 2001) shows that the fraction of 
baryons in the warm/hot intcrgalactic medium 
(that causes most of the absorption) is not sen- 
sitive to the inclusion of additional physics. Note, 
however, that the effects of feedback in these simu- 
lations are much smaller than in the extreme feed- 
back models that assume a uniform preheating of 
the intergalactic gas with energy of order of 1 keV 
per particle in order to explain the observed cor- 
relation of cluster X-ray luminosities and temper- 
atures. Such feedback would affect the gas evolu- 
tion and its spatial distribution lowering gas den- 
sity in groups and poor clusters and completely 
destroying most of the filaments. However, we be- 
lieve that there is no compelling evidence that the 
uniform preheating actually occured. Moreover, 
the existence of lyman alpha forest at high and 
low redshifts argues against such extreme uniform 
feedback. 

The results presented in this paper are only one 
example of useful application of the constrained 
simulations. In Paper I we addressed the prob- 
lem of the coldness of the observed local pecu- 
liar velocity field of galaxies and showed that the 
properties of the local flow can be well reproduced 
in the currently popular ACDM cosmology. The 
constrained simulations are unique because they 
allow properties of individual simulated structures 
to be compared to the properties of their real coun- 
terpart. As a first step in this direction, we are 

5 http: //const ellation.gsfc.nasa.gov/ 



currently carrying out high-resolution constrained 
simulations of the Virgo cluster in the ACDM cos- 
mology. The simulated cluster forms in the en- 
vironment very similar to the real Virgo. The 
differences in the properties and internal struc- 
ture of the simulated and observed clusters should 
thus give us clues about the physical processes 
shaping the evolution of intracluster gas. Another 
very interesting application of the constrained N- 
body+gasdynamics simulations would be predic- 
tions on the spatial distribution of the local ly- 
man alpha forest. These predictions could be com- 
pared to and provide interpretation for the ob- 
served large-scale distribution of the nearby ly- 
man alpha absorbers (e.g. Stocke et al. 2000). We 
conclude that the constrained simulations of the 
type presented here should have many other useful 
applications in comparisons of model predictions 
with various observations of the nearby Universe. 
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